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FRACTIONAL HARDY-LIEB-THIRRING AND RELATED 
INEQUALITIES FOR INTERACTING SYSTEMS 

DOUGLAS LUNDHOLM, PHAN THANH NAM, AND FABIAN PORTMANN 


Abstract. We prove analogues of the Lieb-Thirring and Hardy-Lieb- 
Thirring inequalities for many-body quantum systems with fractional 
kinetic operators and homogeneous interaction potentials, where no anti¬ 
symmetry on the wave functions is assumed. These many-body inequal¬ 
ities imply interesting one-body interpolation inequalities, and we show 
that the corresponding one- and many-body inequalities are actually 
equivalent in certain cases. 


1. Introduction 


The uncertainty principle and the exclusion principle are two of the most 
important concepts of quantum mechanics. In 1975, Lieb and Thirring 
[321133] gave an elegant combination of these principles in a semi-classical 
lower bound on the kinetic energy of fermionic systems. They showed that 
there exists a constant Clt > 0 depending only on the dimension d>\ such 
that the inequality 

( 1 ) 

holds true for every function T G and for all A € N, provided 

that T is normalized and anti-symmetric, namely ||'I'||j;, 2 (RdiV) = 1 and 

T(xi, .,xn) = ...,Xj,...,Xi,.. .,xn), Vi 7 ^ j. 

( 2 ) 

The left hand side of © is the expectation value of the kinetic energy op¬ 
erator for N particles, and for every A-body wave function 'k G 
its one-body density is defined by 
N 

p^(x):='y^ / \'i’{xi,...,Xj-i,x,Xj+i,...,XN)\‘^Y\dxi. 

Note that /g px^/ can be interpreted as the expected number of particles to 

be found on a subset Q C in the probability distribution given by |'kp. 
In particular, = A. 

The Lieb-Thirring inequality can be seen as a many-body generalization 
of the Gagliardo-Nirenberg inequality 

( f |Vu(a:)p(ix^ ( f \u{x)\‘^dx\ > Cgn [ dx, (3) 

\Jr‘^ J / JM.'i 
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for u G Note that for d > 3, the Gagliardo-Nirenberg inequality 

([3]) is a consequence of Sobolev’s inequality 

I|V«IIl2(R'') > Cs\\u\\^2d/(d-2)(^^d-^ (4) 

and the Holder interpolation inequality for L^-spaces. Moreover, Sobolev’s 
inequality can actually be obtained from Hardy’s inequality 

<i>2, (5) 

by a symmetric-decreasing rearrangement argument (see, e.g., [ig Sec. 4]). 

All of the inequalities (I3])-(I1|)-(I5|) are quantitative formulations of the un¬ 
certainty principle. On the other hand, the anti-symmetry which is 
crucial for the Lieb-Thirring inequality ([T]) to hold, corresponds to Pauli’s 
exclusion principle for fermions. In fact, inequality ([T]) fails to apply to the 
product wave function 

4'(xi, X 2 , . . . , Xn) = u{xi)u{x 2 ) ■ ■ ■ u{xn) =■ U®^{xi,X 2 , ■■■, Xn), 

which is a typical state of boson^. In this case Pu®n{x) = iV|u(x)p and we 
only have the weaker inequality 
N \ 

j ^p^^N{xf+‘^/^dx, (6) 

which is, however, equivalent to the Gagliardo-Nirenberg inequality (j3|l. 

The discovery of Lieb and Thirring goes back to the stability of matter 
problem (see m for a pedagogical introduction to this subject). It is often 
straightforward to derive the finiteness of the ground state energy of quan¬ 
tum systems from a formulation of the uncertainty principle such as ([3|), 
(|T|) or dSD. However, the fact that the energy does not diverge faster than 
proportionally to the number of particles — that is, stability in a thermo¬ 
dynamic sense — is much more subtle and for this the exclusion principle 
is crucial. It was Dyson and Lenard who first proved thermody¬ 

namic stability for fermionic Goulomb systems, and their proof is based on 
a local formulation of the exclusion principle, which is a relatively weak con¬ 
sequence of ([2|). Later Lieb and Thirring |32] gave a much shorter proof of 
the stability of matter using their more powerful inequality ([T]). 

Recently, Lundholm and Solovej m realized that the local exclusion prin¬ 
ciple in the original work of Dyson and Lenard 13 [26], when combined with 
local formulations of the uncertainty principle, actually implies the Lieb- 
Thirring inequality ©• From this point of view, they derived Lieb-Thirring 
inequalities for anyons, two-dimensional particles which do not satisfy the 
full anti-symmetry dJj) but still fulfill a fractional exclusion. The same ap¬ 
proach was also employed to prove Lieb-Thirring inequalities for fractional 
statistics particles in one dimension by the same authors [38] , as well as for 
fermions with certain point interactions by Frank and Seiringer m- 

Following the spirit in Lundholm, Portmann and Solovej [36] found 
that Lieb-Thirring type inequalities still hold true for particles without any 
symmetry assumptions — and therefore in particular for bosons — provided 

^In general, bosonic wave functions satisfy @ with a plus instead of a minus sign. 
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that the exclusion principle is replaced by a sufficiently strong repulsive 
interaction between particles. For example, they proved that there exists a 
constant C > 0 depending only on the dimension d > 1 such that for every 
normalized function 'I' G and all G N, 




N 

E 

i=\ 



C [ dx. (7) 

JK3 


The appearance of the inverse-square interaction in ([7]) is natural as it makes 
all terms in the inequality scale in the same way. 


The aims of our paper are threefold. 

• We generalize the Lieb-Thirring inequality ([7|) to the fractional kinetic 

operator (—A)^ for an arbitrary power s > 0, with matching interaction 
\x — The non-local property of (—A)^ for non-integer s makes the 

inequality more involved. Nevertheless, the fermionic analogue of this in¬ 
equality (without the interaction term) has been known for a long time in 
the context of relativistic stability [8]. For the interacting bosonic version 
we will follow the strategy of [36], using local uncertainty and exclusion, 
but we also develop several new tools. In particular, we will introduce a 
new covering lemma which provides an elegant way to combine the local 
uncertainty and exclusion into a single bound. 

• We prove a stronger version of the Lieb-Thirring inequality ([7|) with 

the kinetic operator replaced by (—A)® —Cd^s\x\~‘^^ and with the interaction 
|x — for all 0 < s < d/2. Here Cd^s is the optimal constant in the 

Hardy ineqnality |2T] 

(-A)^-Crf,,|x|-2^ >0. 

Our result can be seen as a bosonic analogne to the Hardy-Lieb-Thirring 
inequality for fermions found by Ekholm, Frank, Lieb and Seiringer ffHEi 

dS]. 

• Just as the Lieb-Thirring inequality ([T|) implies the one-body interpola¬ 
tion inequality dSj), the same will be shown to be true for these generalized 
many-body inequalities. For instance, our bosonic Hardy-Lieb-Thirring in¬ 
equality implies the one-body interpolation inequality 



^d,s\x 



\ l-2s/d 
/ 


u{x)\‘^\u{y)\‘^ 
lx — 



Isjd 


>C [ |n(x)|2(l+2*/‘=')dx, 


for u G and 0 < s < d/2. Moreover, we prove the equivalence 

between the (bosonic) Lieb-Thirring/Hardy-Lieb-Thirring inequalities and 
the corresponding one-body interpolation inequalities when 0 < s < 1. Since 
one-body interpolation ineqnalities have been studied actively for a long 
time, we believe that this eqnivalence could inspire many new directions to 
the many-body theory. 

In the next section our results will be presented in detail and an outline 
of the rest of the paper given. 
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2. Main results 


2.1. Fractional Lieb-Thirring inequality. Our first aim of the present 
paper is to generalize ([7]) to the fractional kinetic operator (—A)® for an 
arbitrary power s > 0, and with a matching interaction lx — The 

operator (—A)® is defined as the multiplication operator |pp® in Eourier 
space, namely 


[{-AYfrip) = m .= 


1 


fix)e-^P-^dx. 


The associated space H‘ 




\2 

Il2(K4) 


(27r)^/2 

is a Hilbert space with norm 

:= (u, (-A)®tt), 


|2 

'H‘ 




and the addition of a positive interaction potential is to be understood as 
the sum of non-negative forms. 

Our first result is the following 


Theorem 1 (Fractional Lieb-Thirring inequality). For all d > 1 and s > 0, 
there exists a constant C > 0 depending only on d and s such that for all 
V G N and for every -normalized function T G 77®(M'^^), 

{"■ (&-"■’■"^ "I 

Since our result holds without restrictions on the symmetry of the wave 
function, and therefore in particular also for bosons, we consider it as a 
bosonic analogue to the fermionic inequalit 50 

^(-Ai)*'h^ > C pv,(x)i+2®/'' dx, (9) 

which holds for wave functions 'k satisfying the anti-symmetry ([ 2 |), where 
the constant C > 0 is independent of N and T. The original motivation for 
such a fermionic fractional Lieb-Thirring inequality has been its usefulness 
in the context of stability of relativistic matter (see [ 8 ] and the recent review 
m ). Our inequality ([ 8 ]) for s = 1/2 and d = 3 is relevant to the physical 


^Throughout our paper, C denotes a generic positive constant. Two C’s in different 
places may refer to two different constants. 
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situation of relativistic particles (which could be identical bosons, or even 
distinguishable) with Coulomb interaction. 

Remark 1. Note that when 2s > d, any wave function in the quadratic form 
domain of the operator on the left hand side of ([ 8 ]) must vanish smoothly 
on the diagonal set 

A := {{xi)fLi G : Xi = xj for some i ^ j}. 

When d = s = 1, it is well known [18] that any symmetric wave function 
vanishing on the diagonal set is equal to an anti-symmetric wave function 
up to multiplication by an appropriate sign function, and hence ([5]) boils 
down to a consequence of Q in this particular case. In higher dimension, 
this correspondence between bosonic and fermionic wave functions is not 
available and it is interesting to ask if a Lieb-Thirring inequality of the 
form Q holds true for all wave functions vanishing on the diagonal set 
(without the anti-symmetry assumption). We refer to Section 13.51 for a 
detailed discussion. 

Remark 2. We have for simplicity fixed the interaction strength in ([5]) to 
unity. One may consider adding a coupling parameter A > 0 to the interac¬ 
tion term and study the inequality 

( 10 ) 

for all A > 2 and all normalized wave functions 'k G with a 

constant C'(A) independent of N and 'k. It is clear that C'(A) > 0 for all 
A, s > 0 and d > 1. However, since the parameter A cannot be removed by 
scaling, it is interesting to ask for the behavior of the optimal constant of 
m in the limits A —5- 0 and A —>• oo. This issue will be thoroughly discussed 
in Section 13.51 

Remark 3. When 0 < s < 1 we can also replace the one-body kinetic oper¬ 
ator (—A)^ by |iV -|-H(x)p® with A G M'^) being a magnetic vector 

potential. By virtue of the diamagnetic inequality (see e.g. m Eq. (2.3)]) 

{u, jfV -I- > (|u|, (-A)®|n|) (11) 

the inequalities l]8 |) -([^- (fT0]l hold with the same constants (independent of 

H). 

When s ^ N, the Lieb-Thirring inequality ([5]) cannot be obtained from a 
straightforward modihcation of the proof of ([7]) in [36] . The non-local prop¬ 
erty of (—A)^ complicates the local uncertainty principle and a fractional 
interpolation inequality on cubes is required. We will follow the strategy in 
[36] . but several technical adjustments are presented. The details are pro¬ 
vided in Section (3] We believe that our presentation here provides a unified 
framework for proving Lieb-Thirring inequalities by means of local formula¬ 
tions of the uncertainty and exclusion principles, and can be used to simplify 
many parts of the previous works [3711381113136]. For comparison, we also 
make a note about fermions and weaker exclusion principles in Section 13.61 
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2.2. Hardy-Lieb-Thirring inequality. Recall that for every 0 < s < fi/2 
we have the Hardy inequalit}!^ [21] 

(-A)* - Cd,s\x\-^^ > 0 on 


where the sharp constant is 

2 , f ri{d + 2s)/4) 
[Ti{d-2s)/4) 


2 


We will prove the following improvement of Theorem [T] when 0 < s < d/2. 


Theorem 2 (Hardy-Lieb-Thirring inequality). For all d > 1 and 0 < s < 
d/2, there exists a constant C > 0 depending only on d and s such that for 
every -normalized) function T E and for all iV E N, we have 



>C [ p^{x)^+^^/^dx. ( 12 ) 

jRd- 

For s = 1/2 and d = 3, the operator in (1121) can be interpreted as the 
Hamiltonian of a system of N equally charged relativistic particles (bosons, 
fermions or distinguishable) moving around a static ‘nucleus’ of opposite 
charge located at x = 0, where all particles interact via Coulomb forces. 

Our result (fT^ can be considered as the interacting bosonic analogue to 
the following Hardy-Lieb-Thirring inequality for fermions: 

(13) 

which holds for every wave function T satisfying the anti-symmetry ([2|) . The 
bound (jl3l) was proved for s = 1 by Ekholm and Frank m, for 0 < s < 1 by 
Frank, Lieb and Seiringer m, and for 0 < s < (i/2 by Frank [15]. In fact, 
m is dually equivalent to a lower bound on the sum of negative eigenvalues 
of the one-body operator (—A)® — Cd^s\x\~‘^^ + V{x) and such a bound was 
proved in [moES]. Unfortunately this duality argument (which has been 
the traditional route to proving Lieb-Thirring inequalities) does not apply 
in our interacting bosonic case. 

Remark 4. The motivation for (|13p was critical stability of relativistic matter 
in the presence of magnetic fields. In both (|12p and (I13p we can, for 0 < s < 
1, replace (—A)® with a magnetic operator |iV -|-H(x)p^; cf. Remark[3j 

The proof of (|13p in [15] is based on the following powerful improvement 
of Hardy’s inequality: For every d > 1 and 0 < t < s < d/2, there exists a 
constant C > 0 depending only on d,s,t such that 

(-A)" - ^ > Cf-\-AY - e on l2(M'^), > 0. (14) 


^The case s > d/2 requires additional boundary conditions at a: = 0 and will not be 
treated here. See m, and m for corresponding fermionic Lieb-Thirring inequalities. 
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Note that by taking the expectation against a fnnction u and optimizing 
over i > 0, we can see that (fH)l is eqnivalent to the interpolation inequality 


uAi-AY- 


Cd,s 
I rrt I 2S 




1—t/s 

'2' >C{u,i-AYu). 


By Sobolev’s embedding (see, e.g., mm for the sharp constant) 

2d 


u, {-AY - 


{u, {—AYu) > Cllrtl'^ 
npl 
Cd,i 


d-2t' 

the bound (fTC]l implies the Gagliardo-Nirenberg type inequality 

\ 


0 <t < d/2, 


X 


2s 


)“> c 


\u 


l—t/s 

2 ' > C||u| 


(15) 

(16) 


2d 


q = 


d — 2t 
(17) 


The bound (flT)) was first proved for s = 1/2, d = 3 by Solovej, Sprensen and 
Spitzer [HI Lemma 11] and was generalized to the full case 0<s<d/2by 
Frank na Theorem 1.2]. 

In fact, ()14p is also a key ingredient of our proof of ()12p . The overall 
strategy is similar to the proof of the fractional Lieb-Thirring inequality 
dS]). However, since the system is not translation invariant anymore, the 
local uncertainty becomes much more involved. We need to introduce a 
partition of unity and use (fTCI) and (HB to control the localization error 
caused by the non-local operator (—A)^. The details will be provided in 
Section m 


2.3. Interpolation inequalities. Let us concentrate again on the case 0 < 
s < d/2. By applying the Lieb-Thirring inequality in Theorem [1] to the 
product wave function T = with ||n||^ 2 md') = 1, we obtain 


N{u,{-AYu) + 


N{N - 1) 


/I 




\x - 

> f dx. 

Jr‘^ 

Since the inequality holds for all iV E N, it then follows that 


^{ 1 *. (-A)“u) + ^ 




> 


<R‘^ \x y\ 

f |^t(a;)|2(l+2s/d)^^ 

jR'i 


(18) 


(19) 


for all /i > 1 (possibly with a smaller constant). On the other hand, by 
using Sobolev’s embedding (fTHI) and Holder’s interpolation inequality for 
L^-spaces, we get 

{u, (-A)*u) > > c YAAA'Y = Y I”!'"*'"''’ P") 

which implies ()19p when 0 < /r < 1. Thus (|19p holds for all // > 0, and 
optimizing over fi gives the interpolation inequality 


{u, (-A)"n)i-2^/'^ 


u{x)Y\u{y)Y 

|x — y| 2 « 


2s/d 
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>C [ ( 21 ) 

for u E ||n ||^2 = 1. Note that in (|21l) the normalization ||u||j ;^2 = 1 

can be dropped by scaling. 

The interpolation inequality ([2111 was first proved for the case s = 1/2, d = 
3 by Bellazzini, Ozawa and Visciglia [1], and was then generalized to the 
general case 0<s<d/2by Bellazzini, Frank and Visciglia [3]. The proofs 
in use fractional calculus on the whole space and are very different 

from our approach using the Lieb-Thirring inequality. 

Remark 5. The inequality (12111 is an end-point case of a series of interpo¬ 
lation inequalities in [3]. The existence of optimizer in this case is open. If 
a minimizer exists, by formally analyzing the Euler-Lagrange equation we 
expect that it belongs to for any e > 0 small, but not 

Thus (I2ip can be interpreted as an energy bound for systems of infinitely 
many particles. 


Remark 6. Note that, when s > d/2, one has 


!L 


R'^xR'i 


|x -y\ 


2s 


dxdy = -|-oo 


for all rt 7 ^ 0 since is not locally integrable. Therefore, the inter¬ 

polation inequality (|21ll is trivial in this case. However, the Lieb-Thirring 
inequality ([5|) is non-trivial for all s > 0 because the wave function T may 
vanish on the diagonal set (see Remark [T]). 


In principle, the implication of a one-body inequality from a many-body 
inequality is not surprising. However, in the following result we show that 
the reverse implication also holds true under certain conditions. 

Theorem 3. For 0 < s < d/2 and s < 1, the Lieb-Thirring inequality ([5]) 
is equivalent to the one-body interpolation inequality m- 

As we explained above, the implication of m from ([8]) works for all 
0 < s < d/2. The implication of dSl) from (I21h is more subtle and we 
obtain it from fractional versions of the Hoffmann-Ostenhof inequality [22j . 
which requires 0 < s < 1, and a generalized version of the Lieb-Oxford 
inequality [nil EH] for homogeneous potentials. We will provide these details 
in Section [5l 


Remark 7. Unfortunately, we can not offer an exact relation between the 
optimal constants in ([8|) and (1211) . On the other hand, from (jlSp it is obvious 
that the optimal constant in ([8|) is not bigger than the optimal constant Ci 
in the inequality 


(/.(-A)*/) + l 


d:,dy >oJ |/(x)|^<‘+^-/-‘) dx. 

k - yr iRd 


for all / E (not necessarily normalized), which is related to the 

optimal constant C in (I2ip by the exact formula 


^ = inf 

C t>o 



^—2sld 



-l+2s/d 
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By the same proof as that of Theorem [3l we also obtain the following 
equivalence for the Hardy-Lieb-Thirring inequality (1121) . 


Theorem 4. Tor 0 < s < d/2 and s < 1, the Hardy-Lieb-Thirring inequality 
(HI is equivalent to the one-body interpolation inequality 






> c 


RdxiRd \x — y\ 
jRd- 


2s/d 


( 22 ) 


The interpolation inequality (|22p seems to be new. Note that the impli¬ 
cation of (I22p from dn holds for all 0 < s < d/2 (by exactly the same 
argument as above), and hence (1221) is also valid in this maximal range. 
There might be some way to prove (1221) directly (as in the proof of (I2ip in 
mis]), but we have not found such a proof yet. 

Finally, we mention that our approach in this paper can be used to prove 
many other interpolation inequalities which do not really come from many- 
body quantum theory. For example, we have 


Theorem 5 (Isoperimetric inequality with non-local term). For any d>2 
and 1/2 < s < d/2 there exists a constant (7 > 0 depending only on d and 
s, such that for all functions u E VF^’^®(M‘^) we have 


\Vu\^^dx 


l-2s/d 


|u(x)P*|u(y)P^ \ 2s/d 

' ^ dxdy ) 


|x - y'/ 




> 


C [ dx. (23) 


This inequality seems to be new and it could be useful in the context of 
isoperimetric inequalities with competing non-local term; see |25l Lemma 7.1], 
[241 Lemma 5.2] and |39l Lemma B.l] for relevant results. The proof of The¬ 
orem [5] will be given in Section [5l 


3. Fractional Lieb-Thirring inequality 

In this section we prove the fractional Lieb-Thirring inequality ([5|). We 
shall follow the overall strategy in |36j , where we localize the interaction and 
kinetic energies into disjoint cubes, but we also introduce several new tools. 


3.1. Local exclusion. The following result is a simplified version of the 
local exclusion principle in m Theorem 2 and Section 4.2]. 

Lemma 6 (Local exclusion). For all d > 1, s > 0, for every normalized 
function 'L E and for an arbitrary collection of disjoint cubes Q’s 

in one has 


E 


l<i<j<N 


\Xj — X, 


12s 


T > > 


E 

Q 


2d"|Q|2V'^ 




(24) 
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Proof. The following argument goes back to Lieb’s work on the indirect 
energy m- Since the interactions between different cubes are positive and 
\x — y\ < for all x,y ^ Q, we have 

1 ^ lQ{Xi)lQ{Xj) 


E 


l<i<j< 7 V 


\Xi — X. 


\2s — 


E E 

Q ^<i<3<N 


\Xj — X, 


12s 




E 

Q 


1 


2d^|g|2V<i 


N 


N 


^lQ(xi) -'^tqixi 


\ 2=1 


2 = 1 


Taking the expectation against T and using the Cauchy-Schwarz inequality 


we obtain the desired estimate. 


□ 


3.2. Local uncertainty. Now we localize the kinetic energy into disjoint 
cubes Q’s. For every s > 0 we can write s = m + a with m G {0,1,2,...} 
and 0 < fj < 1. Then for any one-body function u G we have 

{u,{-Ayu)= [ \p\'^'\u{p)\‘^dp= f \u{p)\‘^dp 

jR'i JRd ^ 2 

^ ^ / \p\‘^'^Ylp'i'"^uip)\‘^dp 


Y, ^{D^uy-AYD^u). 

\oL\=m 


The last sum is taken over multi-indices a = (oi,..., a^) G {0,1, 2,... 
with 

d d d 

|a|=^^ai, a! = J^(Q:i!) and ' 

2=1 2=1 2 = 1 

Here we denoted by p = (pi,P 2 j • • • ,Pd) € and r = (ri,..., r^) G the 
variables in the Fourier space and the configuration space, respectively. 

If s = m, we have 




{u,i-Aru)= Ym \D--\> E 

I I a! jR'* I I “• n JQ 

\Oi\=m l/^l—m (.) 


U\ 


(25) 


\a\=m Q 


for disjoint cubes Q’s. On the other hand, if m < s < m -t- 1, then using the 
quadratic form representatioifl (see, e.g., |14[ Lemma 3.1]) 


{f,i-AYf) = c,,^ [ [ 

jRd. J]Rt 


Ifjx)- fivY 

\x — y|'^+2o- 


dxdy, 


(26) 


‘*Note that this formula only holds for 0 < ct < 1. 
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where 


we have 


Cd,a ■ — 


{U, i-Afu) = Cd,a ^ [ 

II a'- Jm 

\a\=m 


r{{d + 2a)/2) 

\D°‘u{x) — L)“tt(?/)p 




\a\=m 


xR^ 1^ y\ 

|Z)“n(x)-L>"u(y)p 

QxQ 


dxdy 


_ 7/|c^+2cr 


\x - y\ 


dxdy (27) 


for disjoint cubes Q’s. It is convenient to combine ([2^ and (1^71) into a single 
formula 

|2 


where the semi-norm llu 


{u, {-AYu) > 

Q 

2 ^ r2/ 


(28) 




H^iQ) of w £ Li'^iQ) on a cube Q is defined by 
7 JQ if s = m, 

dxdy^ if 0 < a < 


|a|=7n 

„ ml \D‘^u{x) — D°‘u{y)\‘^ 

P jJqxQ 

\a\=m 


\x — yY+‘^^ 

The following estimate plays an essential role in our proof. 

Lemma 7 (Local uncertainty). For every d > 1, s > 0, cube Q C and 
ueL^iQ), one has 


I 1,2 .1 /q C 


Q 


u ■ 


2s/d |Q|2s/<i 


\Ul 


(29) 


for a constant C > 0 independent of Q and u. 

Before proving Lemma [TJ let us clarify a technical point concerning the 
Sobolev space H^{Q) = whose intrinsic norm can be defined by 

(see e.g. [U Section 7.36 and Theorem 7.48]) 

|Q|<m ^ 

Here recall that s = m + a with m E {0,1,2,... } and 0 < u < 1. By 
Poincare’s inequality for H^{Q) (see, e.g., |231 Lemma 2.2]) and the elemen¬ 
tary inequality ja — 6]^ > ^Jaj^ — j6j^ for a, 6 E C, we have 

.. 


Cllu] 


H-iQ) - 


E 


D^u-- 






IQI 


\a.\='in 

From the latter estimate and Sobolev’s embedding, it is straightforward to 
obtain the following equivalence of norms 

|2 ^ 


\u 


Ih'*(q) — ll“llii's(Q) + / kl — Cq\W 


Q 


|2 

IfiTQ)’ 


(30) 


for a constant Cq > 0 depending only on the the cube Q. Now we provide 
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Proof of Lem'ma\^ By translating and dilating, that is, replacing u(x) by 
u(A(x — xq)) for a > 0 and xq G it suffices to consider the unit cube Q = 
[0,1]'^. Then, thanks to (I30p . it remains to prove the fractional Gagliardo- 
Nirenberg inequality 

II^IIh»(q)II'“'IIl 2 (q) > C'||tt||2,9(Q), q = 2 + —, 9 = ^ 2 ^ 

for a constant C > 0 independent of u. Since the (unit) cube Q is regular, 
we may apply the extension theorem to H^{Q) (see [U Theorem 7.41] or 
[481 Theorem 4.2.3]) and obtain for any function u G H^[Q) a function 
U G satisfying 

U\q=U, \\U\\\2l^^d^<C\\u\\L'2(Q), ||lAll^s(Kd) < Cllull^/s(Q), 
where C > 0 depends only on d and s. We will show that 

l|C'll?i.,.-,l|f'llh,V, a C||C/||„,,,„ 9 = 2 + |, 9=^. (32) 

and ([5T]) follows immediately. By Sobolev’s embedding (fTUl) 

4s 2d 

ll^lli7e'>(Rd) ^ C\\U\\L^^d^, q = 2 + — = ^ 
the estimate ()32p follows from the following interpolation inequality 

^ l|f^llil«»(R<^)> V0G (0,1), (33) 

which is in turn a simple consequence of Holder’s inequality 

( [ p‘^"\U{p)\‘^dp\ f f \U{p)\'^dp) > [ p‘^^^\U{p)\^dp. 

/ \JR<i / ./R'* 

□ 

Remark 8. Note that to the semi-norm jj • there is a naturally associ¬ 

ated operator, which for s = 1 coincides with — A^, the Neumann Laplacian 
on H C R'’*. It is a relevant question whether for 0 < s < 1 and bounded 
domains H this operator coincides with (—A^)^ (dehned using the spectral 
theorem), something that was shown in [T3j to be false in the case of the 
Dirichlet Laplacian —A® (see also [laiiaiis] for related results). In any 
case, the analogue of (f2^ for can be proved using the method in 

m- 

We will need the following many-body version of Lemma [3 

Lemma 8 (Many-body version of local uncertainty). For any L'^-normalized 
function T G H^(R'^'^) and for an arbitrary collection of disjoint cubes Q’s, 
the kinetic energy satisfies the estimate 

\ i=l IQ 

where C is the same constant as in Lemma 


1 L c 


C 


(SqP'^) 


.2s/d \Q\2s/d 


pq, 


(34) 
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Proof. Let 7 ^^^ be the one-body density matrix of 'h (see [3U1 Section 3.1.5]), 
which is a non-negative trace class operator on with kernel 

1 ^ f 

'y^^\x,y) := V / ^(xi, ... ,Xj-i,x,Xj+i,.. .,xn)x 

^ JK'i(N-i) 

X .. . ,Xj-i,y,Xj+i,.. .,xn) (35) 

i¥=j 

Since 7 ^^^ is trace class, we can write 

7i^^(3;,y) = '^Un{x)un{y), 

n>l 

where £ L^(M'^) are not necessarily normalized. Then pqt = X]n>i 
and 

N 

i=l 

= (~^) '^n) ^ ) (36) 

n>l n>l Q 

where we have used (|28p in the last estimate. On the other hand, from the 
local uncertainty (I29p we have 


^ = TV [(-A)^7W 




||2 


-V 


c 




d 

d+ 2 s 


^ Q—d/{d+2s) I 


^n\ \\il+2s/d(^Q'j 


for all n > 1. Therefore, by Holder’s inequality (for sums) and the triangle 
inequality we get 



\n>l 


-V 


c 



d 

d+2s 



n>l n>l 

= C <^+^‘>\\p^\\^i+2s/d^Qy 

which is equivalent to 


d 

d~\-2s 


r l+2s/d 

E ll ||2 Jq^'I' 

- C / r 

n>l [JqP'1') 
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The latter estimate and (f36|) imply the desired inequality (IMll . □ 


Remark 9. By using the interpolation inequality (j20l) and the same argument 
of the proof of Lemma [8] (in this case one can work on the whole and no 
partition of cubes is needed), we obtain the following generalization of ([6]): 




l+2s/d 

Pvt 


(37) 


for all normalized functions T E and for a constant C > 0 depend¬ 

ing only on d and s. When 0 < s < 1, (1371) can also be proved using the 
Hoffmann-Ostenhof inequality in Lemma fT^ and Sobolev’s embedding. We 
will use (1371) to obtain the Lieb-Thirring inequality ([5D when N is small. 


3.3. A covering lemma. To combine the local uncertainty and exclusion 
principles, we need a nice choice of the partition of cubes Q’s. The following 
result is inspired by the work of Lundholm and Solovej [371 Theorem 11]. 
In fact, a similar result can be obtained by following their construction. 
However, our construction below is simpler to apply and results in improved 
constants. 


Lemma 9 (Covering lemma). Let Qq he a cube in and let A > 0. Let 
0 < / E L^iQo) satisfy f > A > 0. Then Qq can he divided into disjoint 
sub-cubes Q ’s such that: 

• For all Q, 

//<A. 

JQ 

• For all a > 0 and integer k >2 



• If k = 3, then the center of Qq coincides with the center of exactly 
one sub-cube Q, and the distance from every other sub-cube Q to the 
center of Qo is not smaller than |Q|^/'^/2. 

Note that the simplest choice is A: = 2 and it is sufficient for the proof 
of the Lieb-Thirring inequality dS]). However the case k = 3 will be more 
useful for the proof of the Hardy-Lieb-Thirring inequality (|12l) in Section HI 

Proof. First, we divide Qq into k‘^ disjoint sub-cubes with l/P of the original 
side length. For every sub-cube, if the integral of / over it is less than A, 
then we will not divide it further; otherwise we divide this sub-cube into 
k'^ disjoint smaller cubes with 1/A: of the side length, and then iterate the 
process. Since / is integrable, the procedure must stop after hnitely many 
steps and we obtain a division of Qq into hnitely many sub-cubes Q’s. 

It is obvious that for every sub-cube Q one has Jg / < A and IQI = 

£qp some level i{Q) E {0,1,2,...}. By viewing the sub-cubes 
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Qo 



Figure 1. Example of a division of Qq (in d = 2) with k = 3. 


as the leaves of a full /c'^-ary tree corresponding to the above division (cf. 
Figured]), we can distribute all sub-cubes into disjoint groups {Fi} such 
that in each group F: 

• There are exactly smallest sub-cubes within Fi. 

• The integral of / over the union of these k'^ smallest sub-cubes is 
greater than A. 

• There are at most {k'^ — 1) sub-cubes of every other volume. 

Now we consider each group Fi. Let rrii = inf^gj-. \Q\ denote the minimal 
volume occuring in the group. By the Cauchy-Schwarz inequality we have 



Here in the last inequality we have used the lower bound 

E 

i 5 1QI —2 


/> A> 


>Q 


k^K 

2a 


and that the function t i—>■ t^/k^ — {A/a)t is increasing when t > k^A/{2a). 
On the other hand, using the obvious lower bound 




E 


k^-l 

{k^^rriiF 


A2 k'^-l 

— l)mf 


(40) 
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Here in the second inequality we have used the fact that in Ti , each sub-cube 
has volume k^^rui for some j E {0,1,2,...} and there are at most {k'^ — 1) 
sub-cubes of every volume larger than rrii. Adding (l3^ and (HOl) . we find 
that 



k^-l \ 
4a2(fe'^" - 1 )) 


= 0 , 


where the last identity follows from the choice of a. Since the latter inequal¬ 
ity holds true for every group J-i, the conclusion follows immediately. 

For k = 3 (or any odd integer) there is at each level in the above division 
exactly one cube Q with its center at the center of Qq, and the statement 
follows by iteration. □ 


3.4. Proof of the Lieb-Thirring inequality. Now we are able to give a 
proof of the Lieb-Thirring inequality ([8]). 


Proof of Theorem [IJ By a standard approximation argument we can assume 
that p^, is supported in a finite cube Qq C For every A < p^/ = N^ 

by applying the Covering Lemma [9] with f = p\^, k = 2 and a = 2s/d, we 
can divide Qq into disjoint sub-cubes Q’s such that fq pqr < A for all Q and 


E 


1 

igPd 



> 0 , 


(41) 


with 




Next, from Lemma [U Lemma [8] and (I4ip . it follows that 


N 


2E 

Q 


1 E(-^‘)‘+ E 

r l+2s/d 

ifnP^ C 


Xj — X. 


12 s 




/ 

Jo 


P^ + 


1 


-CA^^/dL 


A 1 


P^ 


IQ 


p^ 




(42) 


for every 0 < A < Ai and for some constant C > 0 depending only on d > 1 
and s > 0. Here in the last inequality in (14^ we have used fqP^ < A for 
all cubes Q’s. 

Finally, using (|42]i for A = {2d^C + l)a =: Aq if > Aq, and using (l371l 
if A^ < Aq, we find that 




N 

E(-'^<)'+ E 

i=l 



l-\-2s/d 


for a constant C > 0 depending only on d and s. The proof is complete. □ 
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Remark 10. Note that, in the case that a coupling parameter A > 0 is 
introduced as in (flOl) . a straightforward adaptation of (H2]) yields C'(A) = C 
for A > 1 and C'(A) ~ for A < 1. 


Remark 11 (Explicit constant). It is possible to derive an explicit constant 
C* in dS]). Let us consider for example the case s = 1 and d = 3. By the 
Hoffmann-Ostenhof inequality (see Lemma fT5]l and Sobolev’s inequality, 

N \ / . N 1/3 r 5/3 


> Cs 




(Irs P'I>) 


2/3 ■ 


Moreover, combining the Hoffmann-Ostenhof inequality and the Poincare- 
Sobolev inequality 

2 


> Cp 


u — 


IQiIq"" 


LHQ) 


as in we get 
/ N \ 


i i=l / Q 

~\Q\ ^ 


2 

LHQ) 


>CpJ/ 

Q 


Q L 


Cp{l-e) 


IQ 


i«(Q) 




'Q 


-2/3 


- Cpie-^ - 1 ) 


|g|2/3 


p^ 


IQ 


for any e G (0,1). From these kinetic lower bounds, following the above 
proof of Theorem [H we find that ([8|) holds true with 

min{(l -e)C'p,C's} 


C = 


Here we can take 


,2/3 


Ao = a(l 6C'p(e - 1)). 


Cs ^(271^)2/3, Cp ig(i^32/3)2(2vr)4/3 


and a = 4 -|- 


VlM 


(the sharp value of Cs can be inferred from [2J HT] and the value of Cp is 
obtained by following m Lemma 1] but it may not be optimal). Then 
optimizing over 0 < e < 1 shows that ([8]) holds true with 


C = 0.002384. 


Although this explicit constant is far from optimal, it is already a significant 
improvement over |36j . 

3.5. Coupling parameter and optimal constant. Let us here consider 
the behavior of the optimal constant of (jlOP as a function of the coupling 
parameter A, 


Celt (A) := inf inf 

N>2 ^ 

I|4'||2 = 1 


r l-|-2s / d 


A > 0, 
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where 'h is in the form domain 

«liV I'P e : f M4{x) := I 

Note that the parameter A cannot be removed by scaling and we are 
interested in the behavior of the optimal constant of (IlOp in the limits A —?■ 0 
and A —oo. We have 


Proposition 10. The optimal constant CbltCA) is monotone increasing 
and concave as a function of X, and satisfies the following: 
i) For all A > 0, any d>l and all s > 0 we have 

0 < min{l, < C'blt(A) < Cgn, 

where > 0 is a constant independent of A and Cgn is the optimal 
constant of the one-body fractional Gagliardo-Nirenberg inequality, 

r ■= inf 

Il“ll2 = l 

a) We have, for all d> 1 and any s > 0, 



lim Cblt(A) = Cblt(O). 

A-J-O 

Moreover, for 2s < d we have Cblt(A) y2s/d 

as A —)• 0, and in 

particular Cblt(O) = 0. 


In addition, we believe the following to be true; 

Conjecture. The optimal constant Cblt(A) also satisfies: 

Hi) Cblt(O) > 0 for 2s > d. 

iv) For all d > 1 and any s > 0 we have 

lim Cblt(A) = Cgn- 

A^oo 

The proof of Proposition [TOl whl be given below. For 2s < d, the limit 
A —)• 0 corresponds to the situation of non-interacting bosons, and by taking 
the trial wave functions T = one can see immediately that Cblt(A) —>■ 
0. However, for 2s > d the situation is more difficult because any wave 
function in FL^i pf must vanish on the diagonal set 

A = G : Xi = Xj for some i j] 

and in particular the trial wave functions are not allowed. 

When d = s = 1, the operator in (fTU|l is that of the Calogero-Sutherland 
model lass], and the limit A —)• 0 on the space Tsym of symmetric wave 
functions is actually equivalent to non-interacting fermions. In fact, Fi) ^ H 
Lgym = \ A) n (see [351 Theorem 2]) and it is well known 

m that any such wave function vanishing on the diagonal set is equal to 
an anti-symmetric wave function up to multiplication by an appropriate 
sign function. Therefore, Cblt(O) is exactly the optimal constant Clt of 
the fermionic Lieb-Thirring inequality ([T]), which is conjectured |33j to be 
Cgn = vrV4. 
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When d = 1 and s = 2, the condition of anti-symmetry is however 
not strong enough to ensure that the wave function is in the quadratic 
form domain which can be seen readily by taking the two-body 

state 'I'(xi,X 2 ) = C{xi — ^ ^i 2 - expect 

C'blt(O) > Clt because of the more restricted domain. 

For d > 2 the situation is yet more difficult: Because of the connected¬ 
ness of the configuration space \ A there is no simple boson-fermion 

correspondence for functions vanishing on A, for any s > 0. Furthermore, if 
s — d/2 G {0,1, 2,then the interaction operator Ws cannot be controlled 
by the kinetic operator by means of the Hardy inequality (see 

[l3l HH]), which makes it difficult to compare T-L^ with \^)- It is 

an interesting open question to determine the complete behavior of C'blt(O) 
in the general case 2s > d. We expect C'blt(O) > 0 for 2s > d because in 
this case \ A) / (by Sobolev embedding), and a smooth 

vanishing condition for 'k on A should imply a non-trivial local exclusion 
principle. In the critical case 2s = d it may happen that C'blt(O) = 0, as 
can be seen for d = 2, s = 1 using the ground state of a gas of hard disks in 
a dilute limit [35] . 

On the other hand, in the limit A —>■ oo of strong interaction, we expect 
the inter-particle distance to go to infinity, and hence the optimal constant 
should tend to the one-body constant Cqn of (|i3]) . It seems that proving this 
would require a concentration-compactness method for many-body systems 
which is not available to us at the moment. We also notice that in the 
physically most interesting case d = 3 and s = 1, the conjectured optimal 
constant in the fermionic Lieb-Thirring inequality ([T|) [33| is strictly smaller 
than Cgn- 


Proof of Provosition UIK We first note that A i—>■ Cblt(A) is the infimum of 
monotone increasing affine functions (denoting T := 

r l+2s/d^^c- l+2s/d’ 

JlRd d>Ir JjRd P^l 

and hence monotone increasing and concave. 

Proof of (i). From Remark HOI we obviously have 

C'blt(A) > min{l, > 0, 

so it remains to prove that C'blt(A) < Cgn- Following [351 Theorem 19], 
we take a sequence of trial states 

'^N,r{x) ■ ^ ^ ^o-(1)(®i) ■ ■ ' £ '^d,N hi ^-‘syvai 

■ itGSn 


with 


Ui{x) := u^{x - Ryi), 


where i?/3)) is a minimizing sequence of L^-normalized func¬ 

tions for (|43p (s.t. both numerator and denominator remain finite), and yi 
are N disjoint points in with \yi — yj\ > 1 for z / j. Since the supports 
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of the tij’s are disjoint, one readily computes that 

, , N (u^J-AYu^) + XCN'^R-'^^ , , 

'^blt(A) < ivJj^,|«R|2(i+2./d) ’ 

and the right hand side of (I44p converges to Cgn in the limit i? —)• oo. Note 
that we could also have taken ^n,r as an anti-symmetric state (a Slater 
determinant). 

Proof of (ii). We will first show that for any d > 0 and all s > 0, 
limA^o C'blt(A) = CbltCO), with 


C'blt(O) := inf inf 


N>2 .r r J+2s/d' 


r i-l- 
jRrf d'F 


(45) 


To do so, we first pick a minimizing sequence ('i^Nk)keN for ([45]1 (with each 
and normalized). Next, we have 


^jVfc,(T +AW,)^jVfc 

0 < C'blt(A) — C'blt(O) < - - — C'blt(O) 


r 

Jr'* P'^ 


l+2s/d 


= ^—r - l+2s/d + —r -" CbLt(O). 


"^k 


(46) 

Given any e > 0, the last term of (|46l) is clearly less than e for k G N suffi¬ 
ciently large, while the first term remains bounded. With such k fixed, we 
then choose A < e(/jjd (^TVfc, so that C'blt(A)-Cblt(0) < 

2e. 

In the case 2s < d we have Cblt(A) ~ as A —>■ 0, which can be seen 
by taking a bosonic trial state 'k = G ^ and letting N x-\ □ 

3.6. A note about fermions and weaker exclusion. In this subsection 
we explain how to adapt our above proof of Theorem [1] to show the fermionic 
inequality ([9]) 

for all d > 1 and s > 0, where the wave function 'k satisfies the anti¬ 
symmetry ([2]) . In this case the kinetic energy not only contributes to a local 
uncertainty principle as in Lemma [8l but also to a local exclusion principle 
of the following weaker form: 

Lemma 11 (Local exclusion for fermions). For any d > 1, s > 0 there is a 
constant C > 0 depending only on d and s such that for all N gN, for every 
L'^-normalized function T G satisfying the anti-symmetry m, and 

for an arbitrary collection of disjoint cubes Q’s in 
IN \ 


>Y. 


c 


\Q\2s/d 


\ i=l / Q 

where q := ^{multi-indices a : 0 < |a| < s}. 


IQ 


p\ii{x) dx — q 


(47) 
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Proof. First, consider one-body functions u E H^{Q) where s = m + a, 
m E N, (T E [0,1). In the case that 0 < ct < 1, we have the fractional 
Poincare inequality (see, e.g., [23l Lemma 2.2]) 


lull ■ > 


^ y 




2 

L^{Q) 


while for jaj = m we have (by iteration of Poincare’s inequality) 

^ ° 0 < |/3| < m. 


Note that JqD°^u = (l,rQ,u) = (r^l,u), where the operator u i-A Ta(u) := 
jaj < m, is relatively bounded w.r.t. the form domain H^{Q). Hence 
we can treat these orthogonality conditions by considering the q-dimensional 
subspace := span{T*l : 0 < jaj < s}. On H^{Q) n we then have 

- |g|2s/dll^lli2(Q)> 


and in general, by taking out the projection onto V^, 

(-A)*|h»(Q) > - ^Vj- 

Now we proceed as in Lemma El although because of the anti-symmetry 
of 'L, the one-body functions Un all have norm less than unity (again, see 
e.g. 1301). We then obtain 


N 




Un 


\H-{Q) - 


2=1 


r2>l Q 


>E 

Q 


C 


\Q\2s/d 


El 

n>l 




which proves the lemma. 


□ 


We note that the Covering Lemma [9] can be also adapted to apply to the 
weaker form of the exclusion principle. This could be useful not only for 
fermions but also in situations when other types of interactions are present 

(cf. 13ZIIIZ1E51ES]). 


Lemma 12 (Covering lemma with weaker exclusion). Let Qq he a cube in 
M'’* and let 0 < f ^ L^{Qq) satisfy Jq^ / > A > 0. Then Qq ean be divided 
into disjoint sub-cubes Q ’s such that 

• For all Q, 

//<A- 

JQ 

• For all a > 0, q >0 and integer k >2, 


y J- 



> 0 , 


(48) 


where 


qk'^\ k'^°‘ — 1 
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• If k = 3, then the center of Qq coincides with exactly one sub-cube 
Q, and the distance from every other sub-cube Q to the center of Qo 
is not smaller than \Q\^/‘^/2. 

Proof. We proceed with the same division procedure as in the proof of 
Lemma m Instead of (p9]) we have 


E 

15 1QI — 


1 



and instead of (|10]) we have 



1 


> - 

mf 


((1 - b)A - qk^^) , 


(49) 






Hence, 


\QY 


f-q 


UQ 


b f]>-bA 


'Q 


QGJi,|C3|>m, 




> - 


j>l 


k'^ - 1 
{k^^miY 


bA k'^ — 1 
mf k‘^°‘ — 1 


• (50) 


E 

Q&Ti 


1 


\QY 


f-q 


UQ 


-b f /) >^{A-qk^-bA{l + 


IQ 




r - 1 
- 1 


from which the lemma follows. 


□ 


From the local uncertainty in Lemma [HI the local exclusion in Lemma 
m and the Covering Lemma \T2[ one can prove the fermionic Lieb-Thirring 
inequality ([9]) by proceeding similarly as in the proof of Theorem [TJ The 
details are left to the reader. 


Remark 12. From Lemma El and the elementary inequality (o^ — a)+ > 
(a—1)+, a > 0, we obtain the following analogue of (|171l for pair-interactions; 


E 

l<i<j<N 


\Xi — X. 


12s 



> 


E 


1 



(51) 


for every normalized function T € lYiW^^). In our proofs of the Lieb- 
Thirring inequality (l8|) and the Hardy-Lieb-Thirring inequality m pre¬ 
sented later, we can certainly use (|5ip instead of (j24p (we then obtain similar 
inequalities but with worse constants). 


4. Hardy-Lieb-Thirring inequality 

In this section we prove Theorem [2j We will need to strengthen the local 
uncertainty principle in Section [3| to account for the Hardy term, and to do 
this we also need a localization method for fractional kinetic energy. 

4.1. Local uncertainty for centered cubes. The following local uncer¬ 
tainty principle is crucial for our proof. 



















FRACTIONAL HARDY-LIEB-THIRRING AND RELATED INEQUALITIES 


23 


Lemma 13 (Local uncertainty for centered cubes). For every cube Q C 
centered at 0 , we have 


|2 



u(x)P 


dx > 


1 /q 


g f I |2 


(52) 


/or a constant C > 0 depending only on d>l and s > 0. 


Note that this local uncertainty principle is significantly stronger than 
the one in Lemma [7] because the left side of (j52p can even be negative. Our 
strategy is to replace u by yu where y is a smooth function supported in 
a neighborhood of the origin, and then apply the Hardy inequality with 
remainder term for yu € To implement the localization procedure, 

we also need the following lemma which controls the error terms. 


Lemma 14 (A fractional IMS localization formula). Let Q be a bounded 
open domain in with d> 1. Let y, 77 : —>• [0,1] be two smooth functions 

such that y(x)^ + r]{x)'^ = 1 and y is supported in a compact subset of H. 
Then for every s > 0, there exists t G [0, s) and a constant C > 0 such that 
for every u G iL®(H), 




Remark 13. It will be clear from the proof of Lemma flTl (provided below) 
that if s G N then t = s — 1, and if s = m + ct with m G {0,1,2,...} and 
0 < fj < 1 then we can take t = s — £ for any 0 < e < min{(T, 1 — it}. 


Note that such a localization bound is well known when 0 < s < 1. In 
the simplest case s = 1, thanks to the IMS formula (cf. [3 Theorem 3.2]) 

|Vu|2 = |V(yu)p + |V( 77 u )|2 - (jVyp + \Vv\^)\u\\ 

we obtain the estimate (|53l) (with t = 0) immediately: 




H^(n) 


^i(Q) llx^^ 

When 0 < s < 1, the estimate 


|2 

'HRn) 


- ||yu 


= [ (|Vy|2 + |Vr7|2)|n|2<C / 
J ^ J 


P- - IItoIP- 


< C 


follows from the representation (I26p 


IlnlP. =Q [[ 

Ix-yl'^+^s 

and the elementary identity (which goes back to a suggestion of Michael 
Loss and was used in [M] l 


|y(x)n(x) - y(y)n(y)|^ + |r 7 (x)u(x) - ri{y)u{y)f - |u(x) - u{y)\'^ 

= [{x{x) - x{y)f + iv{x) - v{y)f] iR[u{x)u{y)]. (54) 

However, the proof of (I53p for s > 1 is rather involved and we defer it to the 
next subsection. In the following, we will give a proof of Lemma [13] using 
Lemma [T3| 
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Proof of Lemma \13i . Since the inequality (|52p that we wish to prove is dila¬ 
tion invariant, we can assume without loss of generality that \Q\ = 1. Let 
X,—)• [0,1] be two smooth functions such that x^(x) -|- rj^{x) = 1, 
x(x) = 1 when |x| < 1/4 and x(x) = 0 when jxj > 1/3. By using 
ry^jup/jxp® < and Lemma [HI we obtain for some t G [0, s) 

2 f 

H‘>{Q) Jq |x|2* 

(55) 




■‘i 

JQ 


12s 


dx > llxitl 


for some constant Ci > 0 depending only on d, s, t (and x)- 

Since x has compact support, x^ can be considered as a function in 
Lf^(R'^). Therefore, by the Gagliardo-Nirenberg type inequality (I17p (there 
taking t = s/(l -|- 2s/d)), 


2 ( llx^^ 


|2 


-c. 


d,s 


\xu\^ 

I rjrt I 2S 


dx I > 


1 / 1x^1'^^+'*/''^ 


c 




2s/d 


(56) 


Moreover, by using the improved Hardy inequality m and the norm- 
equivalence (l3np . we find 




- Cd. 




IymP 

'^ ■dx] \\xu\ 


2s 


1 „ 


|2 

Iff* 


2(l-t/s) 

L2(Kd) 

1 


— ^IIX^IIiyt(Rd) C'llX' 


which by Young’s inequality implies that 


|x|2* > C'llIX'^^lljyqRd) C'||x'w||i2(Rd), (57) 

with Cl as in (f55|) and a (large) constant C > 0 depending only on d, s,t. 
For the function rju, by the local uncertainty in Lemma [71 


2\C^\\h-(q) - Q 


fa 2 

^-cWwWmQY 




(58) 


By using the extension and interpolation arguments as in the proof of 
Lemma [71 we obtain 

\\du\\%\Q)\\w\\\^Q] > C\\vu\\hPq), 
which, together with the norm-equivalence (|3nD . gives the estimate 


2\\vu\\%s^q) > Ci\\7]ufHt(^Q) - <7||r/u||i2(Q) 

for a (large) constant C > 0 depending only on d, s, t. 

By summing inequalities (|5^ - (|56]) - (|57P - (|58P - (|5^ . using 


( 59 ) 
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and estimating the denominators, we arrive at 




,s [ 
JQ 


\x 


„ dx > 

2s — 


c 


Q 


\u\ 


2sld 




for a (large) constant C > 0 depending only on d,s. The final conclusion 
then follows from the elementary inequality 

.2 I „2\P 


+ =2‘ 


= 2‘-P P=l + ^>1. 

d 


□ 


4.2. Proof of the fractional IMS localization formula. 


Proof of Lemma im Step 1. We start with the case s = m E N. Recall 
that in our conventions 

ml 


2 1^. 
\a\=m 


“up. 


Let us consider an arbitrary multi-index a with |a| = m. Using 


D^{xn) = xD^upY. 


a: 


fi<a 






and a similar formula for D°^{riu), we find that 

+ l^"(??w)P = (X^ + ^^)I-D“up 


+ 


E 

f}<a 

+ 2»E 

I3<a 


/3!(a-/3)! 


D^-^xD^i 


+ 


E 

P<a 


(60) 


a\ 


/3!(a-/3)! 


D^-^rjD^i 


a\ 


/?!(«-/3)! 


{xD°‘ ^x + ^ri)D°'uD^u. 


(61) 


Here, for two multi-indices a = (oi,..., ad) and fd = (/3i,..., (dd), the nota¬ 
tion fd < a means jd < a, namely fdj < aj for all 1 < j < d, and jd ^ a. The 
first term of the right side of (|6ip is nothing but \D°‘u\^ since ')^ + rf = 1. 
The next two terms can be bounded using the Cauchy-Schwarz inequality 

2 2 


E 

0<a 


a: 


ld\{a-fd)\ 


D^-^xD^u 




E 

j3<a 


a\ 


fd\{a-fd)\ 


D^-f^r]D^u 


<C^\D^u\^. 

/3<a 


Therefore, by integrating dSH) and using the triangle inequality we get 




+2 EE 

|ci|=m /3 <q: 


ai 


(d\{a-fd)\ 


flxfl 


(x^“"^x + r/T»“-^7/)Zl“uD^u 


(62) 


Now we estimate the last term of (|62p . For every a with |q!| = m, we can 
find 0 < a' < a and 1 < j < d such that L>“ = djD'^'. Note that xdl“~^X + 
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has support in a compact subset of O, so by using integration by 
parts with respect to the j-th coordinate we find that 



+ r]D^-^r])D^uD^ 


u 


= - ^ D^'udj 

Therefore, when |/3| < m - 
estimate 


X + ^r])D^v^ 

X + r]D^~^r])D^u + {xD°‘~^x + 'riD°‘~^vi)djD^v^ . 
- 2, by the Cauchy-Schwarz inequality we can 



X + TjD'^ ^r])D°‘uD^u 


< C\\u 


2 


On the other hand, if /3 < a and |/3| = m — 1 = |q!| — 1, then Z)“ ^ for 
some 1 < /c < d, and hence 

XD^-^X + vD°‘~^v = ^dk (x^ + V^) = 0- 

Summarizing, (1621) can be simplified to 




Since 




Ec 


right side of (|63p by induction and finally arrive at 


— (63) 

, we can continue estimating the 


11 if'" (O) 11 ^'“11 if'" (O) ll^llii'"(Q) 


< C 




This ends the proof when s = m G N. 

Step 2. Now we consider the case when s = m + a with m G N and 
0 < fj < 1. Let us start by considering 




’"^11 oil J JnxQ. 

\a\=m 


\D^{xu){x) - D"(xu)(y)|^ 


dxdy. 


\x — y|'i^+2o- 

We will always denote by a an arbitrary multi-index with |q!| = m. Using 
(fUUl) and the identity |a -|- = |ap -|- 23?((a-|-6)6) — |6p (with complex 

numbers a and 6), we have 

\D°^{xu){x) - L»“(xw)(y)P 

= \x{x)D°^u{x) -x{y)D°‘u{y) 


+ E 

f}<a 


a\ 


/?!(«-/3)! 


^x{x)D^u{x) — Z?" ^x{y)D^u 


= \x{x)D^u{x) - x{y)D'^u{y)\ 
a\ 


IE 

y<a 


f<l(Q-/3)! 


+ 2 »E 


a\ 


fi<a 


/?!(«-/3)! 


D^-^x{x)D^u{x) - D^-^x{y)D^u 
D“(xu)(x) - L>"(xw 
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X ("d" ^x{x)D^u{x) — D" ^x{y)D^'<- 


(64) 

Now we estimate the right side of (|64l) with the help of the Cauchy-Schwarz 
inequality. We have 

D°‘~^x{x)D^u{x) - D^~^x{y)D^u{y) 

= D°‘~^x{x){D^u{x) - D^u{y)) + - D°‘~f^x{y))D^u{y) 

< 2|L>“-^x(x)|2|D^u(x) - D^u{y)\^ 

+ 2|Z)“-^x(x) - D^-^x{y)?\D^u{y)\^ 

< C {\D^u{x) — D^u{y)\^ + \x — y\^\D^u{y)\’^ 
for the second term and 

2{d^{xu){x) - D°‘{xu){y)^ (^D'^~^xix)D^u{x) - D'^~^x{y)D^u 

< \x-y\^^\D'^{xu)ix)-D^ixu){y)\‘^ 

+ k - D°‘~f^x{x)Df^u{x) - D°‘~f^x{y)D^u{y) 

< |x - y\^^\D'^{xu){x) - D"(xM)(y)p 

+ C\x - {^D^u{x) - D^u{y)\'^ + |x - y\^\D^u{y)\'^ 

for the third term. Here we are choosing 0 < e < min{cj, 1 — a}. When 
inserting these estimates into (1641) we find 

||L>“(xu)(x) - H»"(xu)(y)|^ - \x{x)D^u{x) - xiy)D‘^u{y)\‘^\ 

< C\x - y\'^^\D‘^{xu){x) - D°‘{xu){y)\^ 

+ C ^(1 + |x - 2 / 1 “^^) - D^u{y)\^ + \x - y\'^\D^u{y)\'^^ . 

y<a 

Integrating second part of the above inequality against the weight \x — 

y\~^d+2a) 

|ZI“(xu)(x) - i:>“(xu)(y)p - |x(a:)H>"u(x) - x(y)-D“ii(i/)P' 


< 


nxH k “ y\ 

\D^{xu){x) - H>“(xu)(y)p 

fixn 


// 

J Jn-. 


dxdy 


dxdy 


+c’E 

0<a 

< C\\D<^{xu) 


1^; — y|d+2(o--£) 

(1 + |x — y\~‘^^) [\D^u{x) — D^u{y)\‘^ + |x — y\‘^\D^u{y)\'^) 


QxQ 


_ 7/1 d-\-2(7 


\x - y\ 


dxdy 


+ C\\u\ 




where we also estimated difference quotients involving D^u in terms of D'^'u, 
\a'\ = m. Combining the above with a similar inequality for D°‘{r]u), we 
find that 

\\D‘^{xu){x) - D‘^{xu ){y)\‘^ - \x{x)D‘^u{x) - x(2/)E'“tt(2/)P' 

/nxfi 


_ qi\d-\-2o' 


\x - y\ 


■ dxdy 
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+ 


\D'^{r}u){x) — D'^{r}u){y)\^ — \r]{x)D'^u{x) — 

_ ,,|ci+2cr 


nxn k “ y\ 

On the other hand, note that as in ()54D . 

\x{x)D°‘u{x) - x{y)D°‘u{y)\^ + \y{x)D°‘u{x) - y{y)D°‘u{y)\^ 

— \D°‘u{x) — D°‘u{y)'\^ 

= ((x(a;) - x{y)f + ivix) - r?(y))^j5RZ)"«(x)Z?“u(y) 

<C|x-y|2(|Z)“n(x)|2 + |D“n(y)|2). 

Integrating the latter inequality against the weight |x — we get 

\xix)D^u{x) — xiy)D^uiy)\'^ + \r]{x)D°‘u{x) — ri{y)D°‘u{y)\‘^ 


dxdy 

(65) 


f2xO 


|x — y|'^+2o' 
\D°‘u{x) — D°‘u{y)\'^ 

_ ^i\d-\-2a 


dxdy 


dxdy 


< C / \D°u\\ (66) 




inxn \x-y\‘ 

From (l65]) - (|66]) and the triangle inequality, it follows that 

\D°^{xu){x) - D°‘{xu){y)\^ + \D°‘{yu){x) - D°‘{yu){y)\^ 

inxn |x-?/|'^+2<^ 

|F)“m(x) -F)“u(?/)|2 


dxdy 


_^|d+2(T 


nxn \x-y\^ 

2 


-dxdy 
|2 


for all \a\ = m. By taking the sum over all a’s with |a| = m, we get 


11 if* (n) 


< C 


2 I II ||2 I II ||2 


Combining this with the estimate 

|2 ^ x^nu,„,ii2 


which follows from the integer case in Step 1, we can conclude that 


This is the desired inequality. 


< C 




□ 


4.3. Proof of the Hardy-Lieb-Thirring inequality. 


Proof of Theorem [3 By a standard approximation argument we can assume 
that is supported in a finite cube Qq C which centers at 0. Let an 
arbitrary 0 < A < A^. By Lemma 0 with f = k = 3 and a = 2s/d, there 
exists a division of Qo into disjoint sub-cubes Q’s such that Jq p^ < A and 



> 0 , 


(67) 
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with 


b :— — I 1 + I /1 + 


1 - 3 


-d 


^da _ I 


Moreover, for every sub-cube Q we have either that Q centers at 0 or that 
infa;6Q |x| > \Q\^/’^/2. 

Now we claim that there exists a constant Ci > 0 depending only on 
d > 1 and s > 0 such that for every sub-cube Q and for every function 
u G H^{Q) we have the uncertainty relation 

1 /q Cl 




Cd,s [ 

JQ 


\u{x)\^ 


12s 


dx > 


Cl 


(■/q 


2s/d |Q|2s/rf 


[ ( 68 ) 

JQ 


In fact, if Q centers at 0, then (1681) is covered by Lemma [T3l On the other 
hand, if 0 ^ Q, then using \x\ > \Q\^/^/2 we have 

12 o2s 


/ 

JQ 




|2s — 


JQ 


IQ 

and (IU5|) is covered by Lemma [71 Using ([55|1 and arguing in exactly the 
same way as in the proof of Lemma El we obtain the many-body estimate 


N 




2=1 


Q 


1 f r* 

1 JQ Cl 

[!qP^) 


Cl 


I 

JQ 


> 


CiA^^/<i 


l+2s/d 

,P^ 1. 

Q 


^s/d \Q\2s/d 
Cl 




\Q\2s/d 


'Q 


P^- (69) 


Here in the last inequality of (l69l) we have used the bound Jq pq> < A for all 
Q. Combining (|6^ . Lemma [Hand (|67)l . we find that 


N 




\ y2=l 



- CiA^-/'^ i 

l+2s/d , 

P<1, + 

E 

Q 

- i 

l+2s/d , 

P<1, + 

ij 

for all 0 < A < A. 



,Xi — Xj 
l<i<j<N J 


2s 




2d^|g| 


2s 


(^(y J ~ (2d®Cl A 1) J 


2d®Ci-l) 


Q 


2d^\Q\ 


2s 


'Q 


p^ 


(70) 


On the other hand, using the interpolation inequality (I17p with 


2d 


Q = 


d — 2t 


= 2 1 + 


~d 


that is t = 


ds 


d + 2s’ 


and the same argument of the proof of Lemma El we obtain the following 
strengthened version of ([371) : 


( 71 ) 
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for a constant C > 0 depending only on d and s. 

Finally, using dTOll with A = {2d^Ci + 1)6 =: Aq if > Aq, and using 
(HU) if A^ < Aq, we find the desired inequality. □ 

Remark 14. Also in this case it is possible to add a coupling parameter 
A > 0 as in m, and a straightforward adaptation of (I7UI) yields for the 
corresponding constant C'(A) ~ min{l, 

5. Interpolation inequalities 

5.1. Equivalence for the Lieb-Thirring inequality. In this subsection, 
we provide a proof of Theorem [3l i.e. the equivalence of the Lieb-Thirring 
inequality ([ 8 ]) and the one-body interpolation inequality ([ 2 T|) . The implica¬ 
tion of m from dS]) was already explained in Section 12.31 and it holds for 
all 0 < s < (i/2. In the following, we show that the interpolation inequality 
m implies the Lieb-Thirring inequality ([ 8 ]) when 0 < s < d/2 and s < 1. 

We will use the Hoffmann-Ostenhof and Lieb-Oxford inequalities, which 
reduce the kinetic and interaction energies of a many-body state to those of 
its density. 

Lemma 15 (Hoffmann-Ostenhof inequality). For every 0 < s < 1 and every 
normalized function T E one has 

(-A)^Vp^). (72) 

The non-relativistic case s = 1 of dig) was first discovered by M. & T. 
Hoffman-Ostenhof [23 . In fact, (1721) is equivalent to the one-body inequality 
(tt, (—A)®u) > {\u\, (—A)®|u|) (cf. the diamagnetic inequality (fTT]) i and it is 
false when s > 1 . See e.g. |30l Lemma 8.4] for a proof of (1721) and further 
discussions. 

Lemma 16 (Lieb-Oxford inequality for homogeneous potentials). For every 
0 < 7 < (i and for every normalized function T E L^((M'^)^), one has 

{-. ^ ^ 5 // (V3) 

for a constant Clo > 0 depending only on d and 7 . 

The case 7 = 1 and d = 3 of ([73j) was first studied in [271129] . The case 
7 = 1 and d = 2 was proved in [311 Lemma 5.3]. A proof of Lemma fTHl 
following the strategy in m is provided in Appendix lAl 

We are now in a position to complete the proof of equivalence. 

Proof of Theorem\^ We prove that (l 2 T]) implies dSj) when 0 < s < d/2 and 
s < 1. By the Hoffmann-Ostenhof inequality (I72p and the Lieb-Oxford 
inequality (1731) . one has 
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> 


^ J Jw 


' Jw- 


l+2s/d 


k - 

for every e G (0,1]. On the other hand, by using Young’s inequality and the 
interpolation inequality (|2T]1 with u = we obtain 

2s' 


1 - 




«<£k3:M«9 


'xR'* I® 2/1 
p^{x)pq,{y) 

^ xR'* 


\x - y\ 


2s 


\ 2s/d 

dxdy I 


for a constant C > 0 depending only on d and s. Thus 

N 

E(-^*r+ E ,2. 


T, 


T ) > 


, 2=1 


\Xi — 

1<2<J'<A^ ^ 


— Cloe 


l+2s/d 


for all e G (0,1]. As 2s/(i < 1, we can choose e > 0 small enough such that 

C^2s/d _ ^ 

Then the Lieb-Thirring inequality ([8]) follows. □ 

5.2. Isoperimetric inequality with non-local term. In the following we 
show how to use our local approach to Lieb-Thirring inequalities to prove 
the one-body interpolation inequality in Theorem [5j 

Proof of Theorem\^ By a standard approximation argument, we can as¬ 
sume that u is supported in a finite cube Qq C M'^. Let f{x) := |n(x)p^. 
For an arbitrary 0 < A < /, we divide Qq into disjoint sub-cubes Q's 

by applying Covering Lemma [9] with k = 2 and a = 2s/d. Thus we have 
< A for all cubes Q’s and 


^IQI“ Uq' 


f] - 


A 


/ 


'Q 


>0, a := 


1 + 1/1 + 


1 - 2-<^ 


22^ - 1 


• (74) 


Similarly to the proof of Lemma [6l by ignoring the interaction energy 
between different cubes and using \x — y\ < ^/d\Q\^^^ for x,y ^ Q, we have 

2 


fix)fiy) 
k - y\ 


2s — 


E 

Q 


QxQ 


fix)f{y) 

k - 2/1 


2s — 


E 

Q 


1 


\Jq 


f 


On the other hand, by the Sobolev inequality (recall that 1 < 2s < d) 

2sd 


\u 


IiuL 2 »(q) > C'||n||i9((g), q — 


we have 




d-2s 

J jl+2s/d 
2s/d * 


(75) 


(76) 


(/q/) 


Hence 
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\U\^^ = 


> 






— 2 s||„,|| 2 s 


\W^,2s(Q) ^ 


CE 


Q 


Q 


Q 


2sld 


and, combining with (I75]l and (I7il) . 

|w(x)p®|tt(?/)P* 


[ \Vu\‘^^dx + 
jR'i 


\x - y\ 


2s 


dxdy 


> 


> 


Cl 

P^s/d 

Cl 

^2s/d 


in —-U 


f |^|2s(l+2s/d) 


+ 


d^a 


- 1 




Thus, if jgd / > d^a, then we can simply choose A = d^a and conclude that 


|Vu| 2 * + 




\2s{l+2s/d) 


\x — yp^ JiRd 

On the other hand, if f < d^a, then using Sobolev’s inequality 

ll^'“llL2ii(Rd) > C'2||u||j;^2i,d/(d-2s)(lRd), Vu G (77) 

and Holder’s inequality we have 


|Vup® > C2\\f\\ld/(d-2s)(^d) > C 2 


( /r-^ •/') 


> 


^2 


In summary, it always holds that 

|u(x)p®|u(y)p® 


|Vup®dx + 


|x - y\ 


2s 


dxdy > 


"^s/d — (^y^sg^2s/d 

minlCi, 02 } 


(d®a) 2 ®/'^ jRd 


/ 

I 


U 


2s{l+2s/d) 


|2s(l+2s/d) 


(78) 


By proceeding as for the Lieb-Thirring inequality in Section 12.31 that is 
rescaling u i—>• fj,u and optimizing over /i > 0 , we obtain the interpolation 
inequality (f23]l . □ 


Appendix A. Lieb-Oxford inequality for homogeneous 

POTENTIALS 

In this appendix we prove Lemma [T 6 j Note that the argument in the 
original papers [271 [29] uses Newton’s theorem and hence only works with 
the standard Coulomb interaction. The following proof is based on the 
strategy of Lieb, Solovej and Yngvason m Lemma 5.3]. 

Proof of Lemma 1761 We start with the Fefferman-de la Llave representation 

POO P 

L L ’ 

where Bji = B{0,R) is the closed ball in and is a constant de¬ 
pending only on d and 7 (see [T 2 | for Coulomb potential, [28l Theorem 9.8] 
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for homogeneous potentials and m Theorem 1] for more general cases). 
Consequently, 


IL 


p^{x)p^{y) 

\x-y\'r 


dxdy = / fniufdu 


dR 


J^d+'y+l ' 


(79) 


where 

and 

where 


E 


i<i<j<n I** 


Ir ■= P'^ * ISfl 

poo r 

Cd,'r / / gR{u)dt 


-T = 


dR 


^<i+7+l 


(80) 


gR{u)-.= l^, ^ lBR{xi-u)tBa{xj -u)^\ . 

\ / 

Using the Cauchy-Schwarz inequality we find that 

gR{u) = ^ -ix)) ^ -i^)^' 


2=1 


2=1 


= \ fR{u) - \fRiu). 

Combining with the obvious inequality gR{u) > 0 we get 
gR{u) > ^/l(u) - imin{/ij(n),/|(u)}. 

Inserting the latter inequality into (I80|) and using (17^ . we conclude that 

^11 — 

ixRd \x-y\ 

dR 


1 - - —^ 

^ Xi — x,'h' / 2 


P^i-)PM dxdy 


(81) 


^ f [ mm{fR{u)J^{u)}du 
^ Jo JRd. 


J^d+'y+l ‘ 


To estimate the second term of the right side, we introduce the Hardy- 
Littlewood maximal function of 

p*iu) := sup , [ p^ix)dx = \Bi\-^ sup 

R>0 J\^_u\<R R>0 It 


Using /i?(u) < |ili|i?“p*(rt), we find that 

, dR 
(u 


dU fB* 

min{fRiu)jR{n)}j^^< fU 
< Ij* {\B,\R^p*{u)) 


J^d+'j+l 

2 dR 


J^d+'y+l 


Jld+'r+l 

poo 

+ / \B,\R<^p*{u) 

Jr^. 


fR{u) 


dR 


^d+7+l 


dR 


^d+7+l 


I'^ll T^d—'y/ */ \\2 I l-^ll n—7 */ \ 

-- R^ ^{p {u}) + - R^^p [u) 

0 — 7 7 


7 
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for all u and for all i?* > 0. Choosing i?* = {\Bi\p*{u)) we get 


f 


min{/R(«)>/R(«)} 


dR 


< 


|5l|l+7/<i(p*(^))l+7/d 


fld+'Y+l 

for all u E Finally, by the maximal inequality (see, e.g. [45 1 p.58]) 

[ {p*{u))^+"'/‘^du < [ p^{u)^+'r/‘^du, 

jRd. J]^d 

where is a constant depending only on d and 7 , we conclude from (18111 
that 


E 

1 

> - 


1 


\xi — Xjy 




P’q,{x)p^{y) 
R'^xR'^ \x — y\'^ 

This is the desired inequality. 


dcd,-yMd^.y I p 11+7/d f ^+i/d 


dxdy - 

27 (d - 7 ) 


P<$ 


□ 
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